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Abstract 

We compute the four-loop contributions to the /3-function and the anoma- 
lous dimension of the field for the 0(A^)-invariant A^- vector model. These 
results are used to compute the second analytic corrections to the correlation 
length and the general spin-n susceptibility. 



1 Introduction 



Non-linear cr-models have been and are being investigated in theoretical physics for 
a variety of reasons: in condensed matter physics the non-linear A^-vector model 
describes the critical behaviour of systems with an A^-component order parameter [jl|, 
^ ; in elementary particle physics two-dimensional a-models serve as a playground for 
testing ideas which are relevant to four dimensional gauge theories: indeed they are 
asymptotically free and can be studied with a weak-coupling perturbative expansion 

i i Bi- 

The simplest example is the a-model where the fields take values in the sphere 
S^~^ and where the action is invariant under global 0{N) transformations. Besides 
perturbation theory, it can be studied using different techniques. It can be solved in 
the N = CO limit 0, ^ and corrections can be systematically computed |]^, |TD|, |TT|. 
Moreover an exact ^-matrix can be computed [|12], |13[ and, using the thermodynamic 
Bethe ansatz , the exact mass-gap of the theory in the limit (3 ^ oo can been obtained 



14, 0, Iq]. The model has also been studied numerically: extensive simulations have 



been performed for = 3 [y, iV = 4 [|19], |g and = 8 111. The results 
for the correlation length agree with the conventional predictions — including the 
nonperturbative coefficient — to within about 4% for A^ = 3 (at C, ~ 10^), 6% for 
A^ = 4 (for ^ ^ 100) and 1% for N = 8 (for ^ ^ 30). The remaining deviations are 
not much larger than the three-loop correction: for A^ = 3 (resp. 4, 8), at the largest 
/3 where Monte Carlo data are available, the three-loop correction is about 3% (resp. 
2%, 0.5%). For this reason we expect the inclusion of the four-loop term to improve 
sensibly the agreement with the conventional predictions. 

In this paper we compute the /3-function and the anomalous dimension of the field 
up to four- loops, thus extending previous work by Falcioni and Treves [^. From this 
computation we obtain the second analytic coefficients in the perturbative expansion 
of the correlation length ^ and of the vector susceptibility x ^i-nd the third analytic 
correction to the ratio x/^^. Using results obtained in we can also compute 
the second correction to the general spin-ra susceptibility. Some technical details 
concerning the computation are reported in Appendix A. A check of the results is 



provided by the 1/A^ results of [0]: in Appendix B we have checked the correctness 



of the large- A^ limit of our results. We have finally compared our four-loop prediction 
to the available data for the correlation length (a much more detailed comparison for 
A^ = 3, which includes also the susceptibilities, will appear in |2^): we find that the 
discrepancy between theory and experiment at the largest f3 today available is now 
reduced to 2%, 4%, 0.5% for A^ = 3,4,8, the four-loop correction being of order 2%, 
2%, 0.5% in the three cases. The remaining difference should be ascribed to higher- 
loop corrections: for N = 8, using the large- A^ results, we have indeed verified that, 
if all analytic corrections were included up to eight loops, the discrepancy should be 
< 0.1%. 
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2 Four-loop RG Functions 



In this paper we consider the nearest-neighbor lattice A^-vector model in two dimen- 
sions. The fields are unit-length spins cr^. G and the hamiltonian is given by 

-^iatt ^ (^2.1) 
x/i 

The partition function is given by 

Z = [ expi-f^n'"'') Yld^'a.Sial - 1) (2.2) 

As it is well known, the perturbative expansion of this model in two dimensions is 
plagued by infrared divergences. We will not discuss this problem here and we will 
adopt the common technique of adding a magnetic field h to the hamiltonian as an 
infrared regulator. Thus, if the magnetic field points in the first direction we have 
the hamiltonian 

X/I X 

The perturbative expansion is then obtained by considering small fiuctuations around 
the direction of the magnetic field. Thus one sets 



l-7r2,7r,j (2.4) 

and expands the Hamiltonian in powers of tt. 

We will now compute the four-loop /^-function and anomalous dimension of the 
field cr. In principle this can be done through a direct lattice computation. However it 
is much simpler to take advantage of the fact that the four-loop calculation has already 
been done for the continuum theory in dimensional regularization This 
allows us to compute the four-loop contribution by performing a lattice computation 
at three loops. The idea is to compute the finite renormalization constants Zi{j3,fia) 
and Z2{(3, iia) which relate the Green's functions in the MS'-scheme and on the lattice. 

More precisely, define r["j^(pi, . . . /3, /i; 1/a) as the lattice n-point one-particle- 
irreducible correlation function for the 7r-field and rj^(pi, . . . ,p„; /?, h\ fi) its coun- 
terpart in the MS scheme. Then the general results of ||28| imply 

rLt(pi, . . . ,P„; /3, h- 1/a) = Zl''' r||(pi, . . . Zfi/?, Z^Zt^\- /i) (2.5) 

It follows that the /3-functions Wi^fS) and the anomalous dimensions 7(/3) in the two 
schemes are related by 

ly^(zrV) = + (2-6) 

7^^(^r'/5) = 7"^"(/5)- W^"^"(/3)^^ (2.7) 
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In general we expand, on the lattice as well as in the ms scheme, 



W 



scheme 



/32 /33 



^..scheme 
W2 



W 



scheme 



and 



scheme scheme scheme 
schemer n\ _ , Jl 1 , Js 

^ (5 (5^ (3^ 



(2i 



(2.9) 



The coefficients wq, wi and 70 are universal in the sense that they do not depend on 
the renormalization procedure and for this reason we have not added the superscript 
scheme. They are explicitly given by 



Wo 

Wi 

7o 



2% 
N -2 

2n 



(2.10) 
(2.11) 

(2.12) 



All other terms instead are scheme- dependent. In ms they are explicitly given by 



25, 26, 27 



w 



MS 



IN^ -A 
4 (27r)3 



MS ^-2 



w 



{27ry- 



12 



(Ar2 _ 22N + 34) + -C(3)(iV - 3) 



7i 



MS 



A,TQ 3 



1^ 



327r3 
1 



1927r4 



(A^- l)(A^-2) 
{N - 1){N - 2)[4(5 - AT) + 3(3 - A^)C(3)] 



(2.13) 

(2.14) 

(2.15) 
(2.16) 

(2.17) 



where C(3) ^ 1.2020569. On the lattice we have |2 

^2 \ 



Wo 



latt 



and 



„ latt 

7i 



„ latt 
72 



A^-2 
72^ 



N 



TT 



TT 57r^ 



- + — - An'Gi (A^ - 2) + 1 + - - 



24 



«7r 
A^- 1 
(27r)3 



TT 



- — + Att'Gi (A^ - 2) 



llTT^ 

~2A 



(2.18) 

(2.19) 
(2.20) 



where d ^ 0.0461636. 
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We will now compute Wg and ■ We must first of all compute the three- loop 
self-energy on the lattice for the 7r-field. The Feynman graphs are reported in Fig. 
(0). We get 



{N-l)h 



1 



{ 30727r: 



r(3iV-5)(5iV-9)W — 



5127r3 



- 23N + 29 - 27r(3iV - 5) log 



1 



.192(2^ 
:Gi(iV-2) 



17) + — (5A^- 11) 



647r 
log — + A3 



1287r' 



ha? 



32 



+ 



P 



3847r3 
1 



2567r3 
1 

An 



{2N -3){N -l)log' 

lliV^ -41A^ + 36 
^ -(5A^^ -42iV + 



32 



27r(iV - 1) log 



647r2 

i?(2A^-3)(A^-2) 



log 



65) - 

ha^ 
32 



3271 
- i?3 



(77V- 13)- 



f2.211 



where 
A3 = 



7687r 



i3N + 
12 



2N-7 



192 

N 



+ 



24 

3A^ + 2 



2567r2 

It. N 
-Vi+ — 



iN-3) + 



7 



307271' 



:(lliV-21)C(3) 



(2.22) 



5A^- 11 



1 



3847r 
6 , Li 



1287r2 15367r3 



(98A^2- 557iV + 729)C(3) 



-J 



48 



„ + ^(N^ ->iN + 13)V, + ^(N- 2)(N - 3)(2V, - K) 
D 12 lb 



3R+V2- 414 + 2V5 + 2Ve + I6W2) 



(2.23) 



Here Vi, . . . ,Vq, Gi, R, J and W2 are finite lattice integrals. Their definition is 
reported in Appendix A.l and their numerical value in table |1[ 

Most of the graphs can be computed with a limited effort: using the table of 
integrals appearing in Appendix A. 2, all graphs but the last one in Fig. (|I|) have been 
reduced automatically using the symbolic language mathematica. The computation 
of the last graph was much more difficult and involved. The whole computation has 
been done independently by the two authors, and many intermediate results have 
been checked numerically. Some technical details can be found in Appendix A. 

We must also compute the same correlation function in the continuum theory. In 
the MS-scheme we get 



35l^(5iV-9)(3iV-5)log»^ 



2 ^ 5127r3 



h 

UN^ - 23N + 29) log^ — 



Figure 1: Feynman graphs appearing in the computation of the two-point function 
at three-loops. 
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Constants 


Gi 


0.0461636 


R 


0.0148430 


J 


0.1366198 


K 


0.095887 


Li 


0.0029334 


Vi 


0.016961 


V2 


- 0.00114 


V3 


0.07243 


V, 


- 0.0013125 


V, 


0.010063 


V, 


0.017507 


Wi 


- 0.0296860 


W2 


0.00221 



Table 1: Numerical value of the lattice integrals appearing in the (3 and 7 functions. 
Their explicit definitions are reported in Appendix A.l. 



5127r 



/i2 30727r3 



;iliV-21)C(3) 



(2.24) 



^—{2N -3)iN- 1) W ^ + ^(llA^^ - 41iV + 36) log^ A+ 



3847r 



+:TT^(5iV' - 42Ar + 65) log + -^Ri2N -3){N- 2) log - 

/i"' Air yU^ 

96^3 -8){N-2)-^{N- 2)R + ^^(gSAr^ - n7N - 151)C(3) 



2567r3 
1 



From these expressions it is easy to obtain the renormalization constants Zi and Z2. 
We expand 



-4\ 



-4\ 



The terms proportional to (3 ^ and (3 ^ have already been computed in 
three-loop terms we get 



(2.25) 
(2.26) 
For the 



'13 



(iV-2)% 3«V_, (iV-2)2 ,aV 
log -\ 777-^ — (5 + 37r)log 



+ 



647r3 
A^-2 



32 



647r3 



32 



1927r3 



3(A^ - 6) - 187r + (3iV - 20)7r2 - 967r^Gi(A^ - 2) log 



32 



5A^-21 (AT -2)2 3(A^-2) 



192 



647r 



647r2 
7 



1 

1927r3 



(Ar-2)(A^ + 10) 



"i^ AT — fi 1 

■^^(3 - vr)/? - l(iV - 2)(iV - 3)(ir - 2V,) 



+ 



N-2 



^23 



4 

(iV-1) 
1 



^Li - 2Fi - (A^ - 2)^2 - 41/^4 + 2V5 + 2Ve + I6W2 



(2.27) 



+ ^[A^ - 2 + 7r(2A^ - 3)] log' 



647r3 
1 

384^ 
A^-2 



2 oV 
32 



-12(A^ - 2) + (3iV - 17)71^ - 967i^Gi{N - 2) log 



2, ,2 



a /i 



64.3 (1 + 6 



Li iV- 1 



12 



Ar-2 



47r 



(Gi - nV2) 



(2.28) 



Then, using ( p.6| ) and ( |2.7| ), we finally obtain 



A^-2f2iV-7 Ar2-4Ar + 5 5A^-9 73A^-164,, , 

+ + TTT^ + -^TT^ + . C(3) 



27r 



96 



N-2 



271 



(iV-2 + 37r)G 



327r ■ 167r2 
Ar-2 



967r3 



N-2 

(3-7r)i?+^^J-^— Li 



1 



+ (iV-3)(iV-2)(K-2l^3) 

o 

A^- 2 



7^ 



+ 



AT- 1 



2n 



{2Vi + {N- 2)V2 + 414 - 2V5 - 2Vg - IQW2 
5N -21 ^{N- 2)2 ^ 3(A^ - 2) ^ A^ - 2 



(2.29) 



192 



327r 



327r2 



167r3 



A^-2 A^-2 A^-2 

{6N + 37)C(3) [2{N - 2) - Stt] d + — (3 - 7r)i? 



1927r3 
A^-6 



47r 



27r 



4g J-^^Li-l(Ar-3)(Ar-2)(ir-2V3) 



+ 



A^- 2 



4 



[(A^ - 3)^1 + 2(A^ - 2)V2 - 41/4 + 2V5 + 2V6 + I6W2] \ (2.30) 



A check of these results is provided by the 1/A^-results of |TT|. In the large- A^ limit 
we get from the previous expressions 



w 



latt 



1 



ATS 

„.latt 

Is 

ATS 



27r2 
1 

2^ 



^-^ + ^(^^ + 4^2-2^3) 



3271 327r3^^ ^ 27r 16 4 2 8 



+ 0(1/Ar) (2.31) 
+ 0(1/A^) (2.32) 



In Appendix B we have checked that (|2.31|) / ( |2.32|) agree with the predictions of the 
1/A^ expansion. 
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3 Long-Distance Quantities 



We will now use Wg"** and 73*^** to compute the second analytic correction to the 
correlation length ^ and spin susceptibility x- 
Let us begin with ^. In general we have 



exp 



dt 



N-2 



l/{N-2) 



exp 



' 27:13 ' 
N-2 



+ 



1 



Wi 



(3.33) 



The constant is non-perturbative and its value depends on the explicit definition of 
the correlation length. For the isovector exponential correlation length ^y^^\ which 
controls the large-distance behaviour of the two-point function (ctq ■ cr^), an explicit 
expression has been obtained using the thermodynamic Bethe ansatz |T^, |1^, p!6| . 
Explicitly 



l/{N-2) 



r 1 + 



1 



2-5/2 



exp 



TT 



2{N -2) 



(3.34) 



Other possibilities are the second-moment correlation length 



f(2) 



(3.35) 



in the isovector channel, or the analogous quantities in higher-isospin channels. For 
instance, in studies of mixed 0{N)/RP^~^ |^B|, the isospin-2 correlation lengths 
associated with the correlation ((ctq ■ cr^.)^ — 1/N) were introduced: the exponential 
isotensor correlation length ^!^^^^ and the corresponding second-moment correlation 
length ^!^\ Using the fact that in the 0{N) a-model no bound states exist |jl2[, we 
have immediately 



C^(exp) 



1^ 



(3.36) 



For the second- moment correlation lengths no exact expressions exist. However, in 



the large- iV limit we get |29 



c 



f(2) 



C.(exp) I 1 



0.0032 



N 



+ Oil/N^ 



1 - 



1.2031 



N 



0{1/N' 



(3.37) 
(3.38) 



Let us now consider the perturbative corrections to the universal behaviour. The 
first coefficient ai was computed in pT|: 



1 



ai 



N-2 



^ + ^-2.G,) (iV-2) + i-^ 
47r 16 ^ 4 48, 



(3.39) 
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We can now compute 02 to get 



1 



02 



192(A^-2)2 
1 



TT 



327r2(iV - 2) 



[7iV - 16 + 27r(iV - 1)] + J: ^ ™ ~ C(3) 



967r2 N -2 



24(iV - 2) 



-12(A^ - 2)(A^ - 3) - 127r(3iV - 5) - (3iV - ll)7r2 + 2'k'^GI 



T[ N — 6 TT vr 



TT 



+- {2Vi + (iV - 2)^2 + 414 - 2V5 - 2^6 - im2) 



(3.40) 



Analogous expressions can be derived for the isovector susceptibility xv = J2x{^o'^x) ■ 
We get 



2-u)i/iog+7o/too 



exp 



+ 



2wi 



70 



g47r/3/(Af-2) 



' 2nP ' 
N-2 



-(Af+l)/(Af-2) 



(3.41) 



The (non-universal) constant cannot be computed in perturbation theory, and no 
exact expression is available at present. We can evaluate in the large- limit. 
Using the 1/A^ results of |TT| we obtain the following expression: 



C - - 
~ 16 



1 + - (4 + 37c - - 37i^ - 7 log 2) + 0{l/N^) 



^ 0.196 [1 - 4.267/A^ + 0{l/N^)] 
where 7_b ~ 0.5772157 is Euler's constant and 

/r(i/3)r(7/6)\ ^ 



7c = log 



Vr(2/3)r(5/6); 



0.4861007 



(3.42) 



(3.43) 



We will also consider the ratio Xv/^'^ as in this case we can compute an additional 
analytic correction. We write 



'2Txfr 

N-2 



-(7V-1)/(7V~2) 



(3.44) 



where Ry = C^j C| is a non-perturbative universal quantity. 
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The explicit values of 61, ci and C2 are reported in [21, 2^. Explicitly 



Cl 

bi 



1 



{N - 2)2 
1 

27r(A^-2) 



.-i(3iV2 - 23N + 31) + ^ - + (iV - 2)^^! 

96 Stt"^ Stt 



TT 



11 



1 + -{N + 1) + — {N 87r'{N - 2)G 



(3.45) 



Here we will compute 62 and C3. They are given by: 



C3 



1 



{N - 2)2 



-^(3iV - 11)2 _^ ^(7Ar2 _ ^ ^ 1 _ ^^^3^2 _ ^ 
1152 9o 9o 



{2N - 1) + — (3AA2 _ ^3^ ^ ^5^) 



57r 



1 73iV-164 
+ 487r2 N-2 ^^^^ 



6(AA-2) 



6(3iV - 4) - 67i{AN - 5) - (3AA - ll)7r2 + Sn^Gj 



TT N — 6 2tT TT 

-2(3 -tt)R + — J - —Li + -(N -3)(K- 2V3) 

^ ' 12N- 2 3 4^ ^' 

+TT {2Vi + {N -2)V2 + 41/4 - 2F5 - 2V6 - I6IV2) 



N -I 



{N - 2)3 I 384 
N -I 



(3.46) 

^ (lOiV^ - 79iV2 + 186iV - 137) + -^{2N - 3)(3iV2 - 26AA + 37) 



487r= 



1927r 

:(3iV2-18iV2 + 387V-27) 



N-2 



1927r3 



{2N'^ + 57AA - 134)C(3) - —{N - 6) J 



48 



+ [(5AA - 9)7r - 2(2AA - 3)] 
Att 

N -3 



+ {N-1) |i-(3 - tt) i? - ^—^{K - 2Vs) 
iV-7,, 1 



16 



+ ^^V^i + 2(-2V; + + ^6 + 8W2] 



(3.47) 



A check of these expressions is provided by the l/A^-expansion. In Appendix B we 
have verified that these expressions are in agreement with the 1/N results of [|TI| . 
Numerically we have 



1 



bi 
62 



N-2 
1 

{N-2f 
1 , 



-0.0490 - 0.0141Ar) 

(0.0688 - 0.0028AA + 0.0107iV2 - 0.0129AA^ 



AA-2' 
1 

{N-2f 



0.1888 + 0.0626AA) 

(0.1805 - 0.0302AA - 0.0080AA2 _ o.oi08AA^ 



(3.48) 
(3.49) 
(3.50) 
(3.51) 
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Cl 
C3 



N-1 

N-2 
N-1 

{N-2f 

N -I 
(AT -2)3 



0.0908 



-0.0316 - 0.0120iV + 0.0149A^^ 



(3.52) 

(3.53) 



(0.0198 + 0.0005iV - O.OOlliV^ - 0.0102iV3 + 0.0041iV^)(3.54) 



Using the results of |^ we can also compute the second analytic correction to all 
non- derivative dimension-zero operators. A suitable basis is given by 



- " Traces'^ 



(3.55) 



where "Traces" must be such that y^°iv:"n jg completely symmetric and traceless. 
These polynomials are irreducible 0(iV)-tensors of rank n and thus they renormalize 
multiplicatively with no off-diagonal mixing. We define the spin-n susceptibility as 



X 



(n) 



^(y„"--'""(o)F„"- 



(3.56) 



Standard renormalization group arguments give pO 

■ 27r/? 



-[2+n{N+n-2)\/(N-2) 



N-2 



oo ^(n) 



(3.57) 



The non-universal constant C^"-* cannot be estimated in perturbation theory. A gen- 
eral expression is available only in the large- A^ limit. We have 



Q(n) 



(2vr)" 
32 "10 



.0-^0 °ifi(27r)2^ ' \h )t}iPi + ^l 



+ 0(1/A^) (3.58) 



Explicitly C^^^ = it/32, CJ^) = tt^R/4. For n = 2 we also computed the first 1/A^ 
correction to get 



C(2) 
X 



TT 

32 



1 + 



A^ 



-2 - vr + 6 log 



vr 36 , 



C'(2)) + Oil/N') 



0.0982 (1 - 3.171/A^ + 0{1/N^) 



(3.59) 



where C'(2) = —0.9375482. Let us now consider the analytic corrections. In |22 
considered the ratio 



we 



X 



in) 



and computed c^"'* and C2^' . Their expression is 



(n) 



' 2nP ' 
N-2 



-n{N+n-2)/{N-2) 



oo An) 



1 + 1^ 



(3.60) 



„(") 



'-2 



1 n{N + n-2) 



4:71 



N-2 



[TT 



1 / („)x2 n(Ar + n-2) 

2 1^1 ) + Ar32 



(Ar-2)Gi-^-i:(3Ar-14) 
87r 96 



(3.61) 
(3.62) 
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N 








'j^(3—loo'p) 


'j^(4—loop) 


3 


1.85 


89.16(21) 


0.743(2) 


0.781(2) 


0.830(2) 


3 


2.25 


1049(7) 


0.861(6) 


0.897(6) 


0.934(6) 


3 


2.60 


8569(92) 


0.901(10) 


0.934(10) 


0.962(10) 


3 


3.00 


94.6(1.6) -lO^ 


0.930(16) 


0.959(16) 


0.981(17) 


4 


2.50 


34.85(9) 


0.911(2) 


0.930(2) 


0.953(2) 


4 


2.80 


86.07(37) 


0.927(4) 


0.945(4) 


0.964(4) 


8 


5.80 


33.41(8) 


0.985(2) 


0.990(2) 


0.995(2) 



Table 2: Ratio of the Monte Carlo results for ^ and the theoretical predictions at 
two loops, three loops and four loops. = 3 data are taken from |T^, = 4, 8 from 
ITT 



Then we get immediately 

6S"^ = 2a, + cS") (3.63) 



4") = 2a2 + a? + + 2aicr (3.64) 



Numerically 



b^l''^ = ^^^^ (-0.0980 - 0.0283iV + 0.0908 n(n + iV- 2)) (3.65) 



An) 

'J2 



^ 0.1401 - 0.0042A^ + 0.0215A^2 - 0.0257A^^ 



{N-2y 

+ n{n + N - 2)(-0.0363 - 0.0187A^ + 0.0149A^2) 

+0.0041 n2(n + A^- 2)^1 (3.66) 



We want now to compare our four-loop prediction with the available Monte Carlo 
data for the correlation length. We define 

n^-'^^-KP) = (3.67) 

where ^MciP) is the Monte Carlo value of the correlation length Q and ^^"'""^(Z?) is 
the theoretical n-loop prediction given by ( p.33| ) and ( p.34| ). In table (0) we report 71 



for some selected values of N and /?. It is evident that the inclusion of the four-loop 
correction improves the agreement between theory and "experiment" . The remaining 
discrepancy at the largest /5- values available today is now 2% for 0(3), 4% for 0(4) 



^ We consider here the isovector exponential correlation length ^^^P' . Notice that the data in ||n 

c(2) 



for = 3 refer instead to The two quantities differ however by less than 0.1% Is^ and thus 



we will ignore the difference. 
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and 0.5% for 0(8) and it can presumably be ascribed to the neglected higher-loop 
corrections. 

We want now to try to keep into account the higher-loop corrections using the 



large- iV results of ||TT|. Let us first define 

aW^'^N) = N^-Hn (3.68) 

The coefficients a^^''^-* are the leading contribution to a„ in the limit N oo and a„ 
are numerical coefficients which can be computed using the expansion. Their 
explicit value for n = 1 ... 6 is reported in Table (|^) . For n = 1 and n = 2 we 
can compare a^^^^ with the exact value a„. One finds that in both cases the ratio 
ttn/a^^^^ is a decreasing function of N reaching the limiting value of one for — > oo. 
The convergence is however slow. For = 4 we have indeed 



while for = 8 we get 



«i(4) 
aS^/^)(4) 

^2(4) 

ar/^^(4) 

ai(8) 
af/^)(8) 

«2(8) 

4^/^) (8) 



3.73 (3.69) 

2.88 (3.70) 

1.91 (3.71) 

1.58 (3.72) 



The approximation is good at the 10% level only for > 50 (resp. A^ ^ 35) for ai 
(resp. 02). Nonetheless we can try to use the 1/A^ results to get a rough idea of the 
role of higher loop corrections. 

A^ = 8 should be the case where the approximation works better. In this case we 
will assume that the coefficients an{N = 8) for n > 3 are given by 

an{N = 8) = fc„8"-W (3.73) 

where a„ are defined in Table (|^) and kn is a number that we will vary between about 
1 and 2. In this way we can get an estimate for n^^-i°°p) , For (3 = 5.80 (see Table (Ef)) 
we obtain for kn = 1 7^(8-'°°^') = 0.998; while for kn = 2 we have 7^(^-'°°P) = 1.001. 
The eight-loop correction is of order 3^6 ■ 10~^. Thus at this order we would expect 
an agreement at the order of 0.1% and this is indeed what we get from this rough 
approximation. 

We can try the same rough approximation for A^ = 4 , writing, for n > 3, an{N = 

4) = A;„4"-ia„. For (3 = 2.80 we get n^^-^°°P^ = 0.970 (resp. 0.994) for fc„ = 1 (resp. 

5) . Although in this case it is very difficult to make any quantitative statement this 
calculation shows that the numbers are in the correct ball-park. 

From this analysis we can thus conclude that the theoretical prediction of [|H|, is 



in very good agreement with the Monte Carlo data. We do not discuss here other long- 
distance quantities like the vector and tensor susceptibility. A detailed comparison 
with Monte Carlo data will appear elsewhere . 
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A Technical Details 



A.l Notations 

In this Appendix we introduce the notations we have used in the explicit computation 
of ^3 and 73. The one-loop perturbative results will be written in terms of the integral 



/TT 
-IT 

where = 2sin(p^/2), = T,^pI and 

A{p) = 

p^ + h 

The integral I{h) is logarithmically divergent for h ^ 0. Explicitly 

2 , / 4 



I{h) 



-K 



7r(4 + hy' \A + h, 
where K is & complete elliptic integral of the first kind. We will also use 



h{h) 



^ (fp 



A{pf 



dl{h) 



(A.74) 



(A.75) 



(A.76) 



(A.77) 



(27r)2 dh 

We need also some basic two-loop and three-loop integrals. To simplify the notation 
let us first introduce 

dPq d^r d^s 



dp{p, h) 



dfj,{h) 



(27r)2 (27r)2 (27r)2 

d'^q d'^r d^s dH 



{27ryS{p + q + r + s)A(g)A(r)A(s) 



(A. 78) 



i2nf 6{q + r + s + t)A(g)A(r)A(s)A(t) (A.79) 



(27r)2 (27r)2 (27r)2 (27r)2 
We will use the following integrals which already appear in our previous work 



R 
J 
K 



lim 



h j du(0, h) 
j dii{Q) [{qTr) 

J dm ^ 



2 -2 -2 
q — r 



[sTt) 



2 52 - £2 



g2(s 



{s + tf -{s + t)^{s + t),f' 



(A.80) 
(A.81) 
(A.82) 
(A.83) 

(A.84) 
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We have moreover introduced a set of 8 new constants . The quantities Vi, . . . Vg 
correspond to lattice infrared-finite integrals and are explicitly given by 



V2 



1 



(A.85) 



d/i(0) 



(q'+'r) 



[s + tf - - ? 



X 



E 



^3 



rf/i(0) 

rf/i(0) 
d/i(0) 



(g'+~r)2 
1 
1 

{q'+^r) 



[sTtf - - Pi [(gTr)^ - g2 - f2 



E ^J^l E sin sin 

H u 

E ^ sin2(g + r)^ 



+ ^j^^M^M COS 



q + r 



(A.86) 

(A.87) 
(A.88) 

(A.89) 

(A.90) 



We introduce also Wi and W2 as the finite part of infrared-divergent integrals: 



W2 



lim 



lim 



J du{0, h)A{q) Y: {qt - 2gjrj) + hihf - ^I{h) 
J dfi{h) A{q + r)^^ sin sin r^^ sin 5^ sin t^. 



(A.91) 



6 



1 1 



57r 



1 1 

+ 



1 



R 



32 167r2 167r 



/(/.) 



(A.92) 



The numerical value of all the constants is reported in Table |T|. 

The numerical computation of Vi, . . ., Vq does not present any difficulty as the 
integrals are infrared-finite. More tricky is handling the integrals leading to Wi and 
W2- In this case we have used a method which was introduced in the context of the 
1/N expansion in [ 3^ . Let us consider the case of Wi . First of all let us determine the 
divergent terms which must be subtracted from the original integral. If we introduce 



di,iq,h) = 

we can rewrite 

[du{0,h)A{q)Y.{ql-2qlfl 



(27r)2 



A(r) A(g + r) (gTr)^f^ 



<fq 



(A.93) 



(2^)2^(^)'ECcos|rf,(g,/.) (A.94) 
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Then we determine the behaviour of d^lq, h) for g — by writing an integral repre- 
sentation for df^{q, h). Using the technique presented in we get 



ZTT 2 Jo af^ [ [aia2)^'^ 

0^(4 + h)-{A + hf + 4x(l - x){2ql - f 

where E{() and K{() are complete elliptic integrals and 

>2 4aia2 



(A.95) 



(aia2)3/2 1 - (2 



(4 + /i)2 



ai — a2) 

2 



1 — x(l — x)q 



Using (|A.95| ) we can now compute the expansion of 
q^/h arbitrary. We get 



h) for q"^ 0, h 



(A.96) 

(A.97) 
0, with 



Sr\q,h) 



h 



Ah" 



— n - log — + - ( 1 + — ) log 



H 

q' 



32 

2h , 
1 + — log 

q p 



where 



Then we rewrite ([A.94|) as 



(A.98) 



(A.99) 



-2 



-2 



d^q 
(27r)2 

d'q 



(q' + hr , 
1 



E«cos|rf,(g,/.)-^^^,^ 



,2 T.id'^r\l^h) 



(A. 100) 



(27r)2 (g2 + h) 

The first integral is infrared finite. Thus we can take the limit h ^ obtaining 



-2 



d^q 
{2ny 



1 



(^2)2E^cos|^d^(g,0) 



lfl-2^-llog^^ 



(A.lOl) 



Although this integral is finite its numerical evaluation is complex as the integrand 
is a difference of two divergent quantities. To get stable results we have split the 
integration domain in two parts: a disk of radius e = 0.1 around the origin and 



the remaining region 



-TT, n] — D^. The integration over is done numerically; 



to compute the integral over we have first expanded the integrand up to 0{q^^) and 
then we have performed the integration analytically. In the implementation a useful 
check is provided by the expression of (i^(g, 0) along the diagonal, i.e. for q = (Z, I) 
which can be computed exactly 



d,{{l,l),0)= -—log 



tan ■ 



(A. 102) 
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Let us now compute the second integral in ( A.lOOp which is still infrared divergent. 
We first change the integration domain: if A > i{\pl we rewrite it as 



-2 



/da-[-^,^]2 (27r)2 (g2) 
The second integral is easily computed numerically, while the first gives 



(A.103) 



1 2 ^ 3 h 

log - log — 



1 



R 



327r 



A2 327r2 °A2 327r2 



In the calculation the choice of A is completely arbitrary. We have chosen A 
as in this case it is simple to replace log/i with lih). 

The calculation of W2 is completely analogous. Introducing 



(A. 104) 

32 



(2^ 



A(r) A(g + r) sinr^ sin(g + r)y 



we must compute 



d^ 



First of all we compute an integral representation for d^y{q, h). We get 

dii{q, h) = cos y di{q, h) - ^ di{q, h) 

where di{q, h) is given in ( [A.95|) and 

4 



(A.105) 



(A. 106) 



(A.107) 



di{q, h) 



2 TT 

— f dx 



A + h 



2cos2 — ((4 + /i)^ + a?-aJ 



- cos' ^(4 + /i)(4 + a\) - (1 - 2x)'(4 + h) 



■a2 e 



2 /-i 



TT Jo 



dxK{Q- 



aia2Y^'^ ai 



-4cos2|- + l(l-2x)'(4 + /i)g2 
2 a2 



(A.108) 



and 



(ii2(g, h) 



A + h 



TT 



■ sm qi sm g2 x 



dxxil — x) 



c 

[aia2)5/2 



2-c: 

1-C 



i5(C) - 2ir(C) 



(A. 109) 



18 



From these integral representations we thus get in the hmit ^ 0, h ^ with q^/h 
arbitrary 



2 

1 QtiQi 



1 + — logK— - 



471 g2 

The computation is analogous to the previous case. The relevant integral is now 



(A.llO) 



r (Pq I ^ .g^p) 



log 



47r V32 167r2 IGtt 



1287r2 V tt; ° A2 
1 i?\ , /i 3 , 

77 log + IT 1 - - + 



A2 8 



TT 
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5767r3 



a'i) + o{h\oih) 

(A.lll) 



where is a disk of radius A around the origin. 

In order to check the manipulations of the extremely cumbersome expression for 
di{q,h) we found very useful the expression oi di{q,h) for g = (/, /) and h = given 

by 



7 //, ,^ ^^ 1 1 1 Sin^m, 



tan ■ 



(A.112) 



A. 2 Some Lattice Integrals 

In this section we report a list of integrals we have used in our computation. 
Two-loop integrals: 



I diy{p,h)A{q)A{p + q)-^ = 
j dp{p,h) A{q)A{p + r)~'^ = 



(A.113) 
+ 0(h',hp',p') 



uh)-l[i{h)-\)+i(hr-\m + 



h 



p 



An 



+ 0{h,p' 



(A.114) 



j du{p, h) A{q) (2 A{p + q)~^A{p + r)"^ + A{p + r)~^A{p + s)"^^ 

h{h) + 21{hf + ^p' [i- [l{h) - i-) + ^ + 0{h)^ + 0{hy) (A.115) 
j du{p, h) A{p + q)-^ = 
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du{p, h) A{p + q)-^ = 



IQtt ^ ^ 12 6 . 



2m-\ + h[-i{hf + \i{h)) 



+p' 



I{hf + I{h) {---]- — + 2Gi-2R 



TT 
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j dv{p, h) A(p + q) ^/\{p + r 
I{h)+p' 



-1 



ZTT 4o 



du{0,h)A{q)J2sm^q^ = 

I J du{0) - hihf - + l(2i? - W,) + 0{h) 



Three-loop integrals: 



7 



^.WA(. + r)-^l/(/.)3-^- ^^^^^ 



C(3)+0(/i) 



y dfi{h) A{q + r)-2 = 4/(/i)2 - /(/i) + K + 



-I{hY + 2(1- -)l{h) 



TC 



J 1_ 77 

^8^ 4^ ^192^^^ ^ 

1 (i/x(/i)A(g + r)-iA(g + s)-i 

2W-l-^ + ^ + 
^ ^ 24 2 6 
r 4 _ 1 



2 2V| 

-Gi + + 2^2 

TT 3 



--i{hy + -i{hy+[AGi- — 



1 



TT 

1 35 



1 



1 



12 27r^ 



/(/.) 



Ta ^ ~ mo s^(^) ^ - ^2 - -Li 

Ibn 8tt'^ 1927T'^ n 3 



A. 3 Analytic evaluation of one-loop integrals 



(A.116) 



0{h\hp',p^) (A.117) 



(A.118) 



(A.119) 



(A. 120) 



0{h^) (A.121) 



(A.122) 



In this appendix we discuss the evaluation of the most general one-loop integral. Let 
us introduce the notation 



B{r;nx,ny 



- d^k 



(A.123) 



for rij > 0, r > 0, and x, y all different. In the following when one of the arguments 
Hi is zero it will be omitted as argument of B. 

Using a technique we developed for four- dimensional integrals ||3^ all these inte- 
grals can be reduced to a sum of I{h) and hih). 
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We will firstly generalize ( |A.123| ) by considering the following more general inte- 
grals 

Bs(r;n^,ny) = - — ^ A.124 

where r is a positive or negative integer and 6 a real number which is introduced in 
order to avoid singular cases at intermediate stages of the computation and which 
will be set to zero at the end. 

The first thing we want to show is that each integral Bs{r; n^, Uy) can be reduced 
through purely algebraic manipulations to a sum of integrals of the same type with 
n^ = ny = 0. 

Indeed the integrals Bs satisfy the following recursion relations: 

Bs{r-1) = ^[Mr-l)-hBs{r)] 
Bs{r;m,l) = Bs{r - l]m) - B5{r;m + 1) - hBs{r]m) (A.125) 
which can be obtained by the insertion of the trivial identity 

P = ^fc; (A. 126) 

5=1 

and by keeping into account the cases in which the index q equals one of the other 
indices. Furthermore, when m > 1 we can write 

,fer fer- (A,127) 

Then, integrating by parts, we obtain the recursion relation: 

Bs{n-...,m) (A.128) 

Tfl — 1 ~ AiTfi — 6 ~ 

-Bsin — 1; . . . , m — 1) Bsiji — 1; . . . , m — 2) 



(n + 5-1) ' ' ' n + 5 

+ ABs{n- . . . ,m - I) (A.129) 

These relations allow to reduce every integral Bi{r\ rij., Uy) to a sum of the form 

r-l 

Bs{r-n,,ny)= J2 Mp)(^piS) (A.130) 

p=T-nx-ny 

For p > 1, lim^^o cip{6) is finite while for p < 1 ap{6) may behave as 1/6 when 6 
goes to zero, meaning that we need to compute Bs{p) including terms of order 6 when 
P<0. 

Now let us show that all Bs{p) can be expressed in terms of ^3^(1) and Bs{2). 
Indeed let us consider the trivial identity 

Bs{r-A,l) - 2Bs{r + l;2,l) - hBs{r + = (A.131) 
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If we apply the reduction procedure to the integrals appearing in this relation we get 
the following two relations 



A + h 



q{3q^ + 3g + 1) - (3g + 2)5 Bsiq + 1) 



+ 



q + l 
q3 



(32 + 2Ah + 3h^){-qiq + 1) + 25) Bsiq + 2) 



(A.132) 



^ h{4 + h){8 + h) ^^^^ + + 2) - 5i3q + 4)] Bsiq + 3) + 0{5') 



Bsiq) 



1 



iq-2)iq-l)hiA + h)i8 + h) 



Bsiq -3) 



3g2 - 15g + 19 



+ 



(g- l)(g-2) hi8 + h) 
q-23h^ + 2Ah + 32 ~ 



Bsiq - 2) 



q-lhiA + h)i8 + h) 



Bsiq-l) + 0i5) 



(A. 133) 



If we now apply the first relation for g < and the second for q > 2 we express every 
integral as 



Bsir; n^.,, n„ 



Ai5)BsiQ) + Bi5)Bsil) + Ci5)Bsi2) + 0(<5) 



(A. 134) 



A careful analysis of the structure of the recursion ( |A.132| ) shows that -8(5) and C(5) 
are finite for 5 — > 0. As the l.h.s. is also finite for 5 — >■ we get that also Ai5) is finite 
in this limit. Thus we can set 5 = to get 



B(r;n,,n,) = A(0) + 5(0)/(/i) + Ci^)hih) 



(A.135) 



As a final remark, notice that the whole procedure is completely algebraic and can 
be easily implemented on a computer using a symbolic language. 



B Comparison with the expansion 

In this appendix we want to compare our four- loop result with the 1/A^ results of 



TT| . In the perturbative limit we have 

X 



eo + ^6 + Oil/N') 
Xo + ^Xi + Oil/N') 



with 



Xo 



1 


27T 13^ 


^^exp 
V32 


[ N J 


N 1 




32/3 


) 
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(B.136) 
(B.137) 

(B.138) 
(B.139) 



and 



11 



47r/5 , 47r/5 
log- 



TT 



Xi 



Xo 



N 
-2vreo 

-2 



N 



21og2-7s + l-- 



lat {27t)^ (3/N + Ao{k) 
d?k 1 



{A^{k)-B^{k)Ao{k)) 
1 



-47rxo 



fc2<32 (27r)2 F log(A;V32) + 47r/3/Ar 

2 - TT 



(B.140) 



47r/? 

3 log — 4 log 2 - 37i5 + 37c 



cPk 



1 



/at {2nY P/N + A^^k) 
cPk 1 



A(A:) - B,{k)A,{k) - ^B,{k) - -i- 



1 



where 



Ao{k) 
A,{k) 



fc2<32 (27r)2 A;2 log(A;2/32) + A7r(3/N 



(fq k"^ -f - {k + qf 



(B.141) 



(B.142) 



and 



(27r)2 (27r)2 
-k\f^ 



B,{k) 



{2Txf5{k-r 



-fii(A;)fV(f2 + s2)}B.143) 



(B.144) 



A point must be noted in the solution (|B.140D , (|B.141D : both and xi ^^e expressed 
in terms of differences of two integrals which have a non-integrable singularity for 
= 0. This notation is however symbolic and it must be interpreted in the following 
way: for two functions f{k) and g{k) we define 



(Pk 

lat 



f{k) 



(Pk 
fe2<32 (27r)2 



9{k) 



dPk 

B (2^ 



{f{k)-9{k)) 



(Pk 
c (2^ 



g{k) (B.145) 



where B is the domain [— tt, Tr]^, C = — B and the disk centered in the origin 
of radius r. It is easily checked that with this definition everything is well-defined for 
= as in both cases the singularity cancels in the difference. From ( |B.140|) and 
( [B.141|) we can easily derive the large- iV behaviour of 02, 62 and C3. We get 



-27T 



+ 



(Pk 
lat (27r)2 
1 



A,{k) {A^{k)B^{k) - A^{k)) 



<fk 1 fc2 

87r2 yfc2<32 (27r)2fc2 



+ 0{l/N) 



(B.146) 
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N 



-An 



"^'^ Ao{k) (Ao{k)B,{k) - A,{k) + ^B,{k) + ^ 



+ 



lat (27r)2 

3 



Airk^ 



(fk 1 , 



+ 0(1/Ar) 



iV2 



C3 



167r2 A2<32 (27r)2 P ^ 32 
Zat {2,7:) \ k^ 

IdTT^ Jfc2<32 ZTT r 32 



(B.147) 



(B.148) 



In order to compare these expressions with our results we must reexpress these inte- 
grals in terms of our basic constants. In the Appendix F of |^ one can find the 
basic relations which are needed to compare the results of the expansion with 
the standard three- loop result. As we will use them in the following we report them 
here: 



d^k 1 



£k 1 



lat (27r)2 k2 yfc2<32 (27r)2 P 

d'^k Ao{k) r d'^k 

lat (27r)2 k2 

d^k 





1 , fc2 







/ 

Jk 



lat (27r)2 

d^k 



lat (27r)^ 



A,{k) + 
Biik)Aoik) 



fe2<32 (27r)2 47rP ^32 

_d^J_(. k^_ ' 
fc2<32 (27r)2 27rA;2 1^ 32 " \ 



87r2 



d'^k 



log 



fe2<32 (27r)2 27rA;2 ^32 



1 

32 



(B.149) 
(B.150) 
(B.151) 
[B.152) 



We must now derive analogous expressions for the integrals appearing in 02, 62 and 
Let us begin by considering the quantity 



d'^k Ao{kf 



1 



d^k 1 



hat (27r)2 A;2 IGn^ Jk^<32 {2tc)^ P 
In this case we start from the three-loop integral 

T = [ rf/i(/i)A(g + r)-i 



log2- 
^ 32 



(B.153) 



(B.154) 



whose value is reported in ([A.120| ). It can also be computed using the method we 
presented in Appendix A.l when we considered Wi and W2- If we define 



we can rewrite T as 



doiq, h) 



T 



d'^r 



■A(r)A(r + q) 



d\ 
(27r)2 



A{q)-'do{q,hy 



(B.155) 



(B.156) 



''Notice the change in notation: ^^-^ = 2A#, B'^' = 2B^ 
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Now for g ^ 0, /i — s> 0, /h fixed, dQ^q, h) becomes its continuum counterpart 



where p is defined in ( [A. 991 ). Thus we rewrite 



4 



A(A:)-ido(fc, hf - / -^(^2 + h)dr{k, hf 

lat [2t:Y Jfc2<32 (ZTT)"' 



The second integral is easily done and gives 

/ + h)dl''^{k, hf = \og^ — ^C(3) (B.159) 

For the term in curly brackets let us notice that 

A{k)-^do{k, hf - + h)dl^%k, hf (B.160) 

is integrable for ^ for all values of h. Thus we can expand the integrand for 
/i — >^ 0. Using the fact that for /i ^ at A; fixed we have 

do{k,h) = ^{I{h) + Ao{k))+ 0{h\ogh) (B.161) 
k^ 

and the relations (p.l49|) /( |B.150| ), we get at the end 

^4|/ ^Wh;] (B.162, 



hat (27r)2 p 167r2 7fc2<32 (27r)2 F ° 32 
Comparing with ( [A.12(j| ) we get 



liat (2Tif A;2 167r2 yfc2<32 (27r)2 P ^ 32 967r3^' ' 12 
Simple algebraic manipulations give also 

d'^k ,,,2 1 /■ ci^A; 1 1 ..on K Vi 



B^ik)Aoikf+— log2 =_ ((3)-- + -i + -i 

lat (27r)2 87r2 7a;2<32 (27r)2 k^ 32 487r'^ 16 6 8 

(B.164) 

The last integral which is needed is 

-Ai(A:)Ao(A;)+ — / log— -1 log— (B.165) 



lati2TTf ' ' 87r2 A2<32 (27r)2 P \ ^32 / ^32 
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This quantity can be handled exactly in the same way starting from 



J dfi{h)Aiq + rr'Aiq) = -2 J ^^A{q)~'do{q, h) 



dh 



[B.166) 



which can be computed by derivation with respect to h from ( A-.121| ). The relevant 
integral is 



(?2<32 (27r) 



jA(g) dQ^[q,h) 



dh 



2TT^h 



log 



h 

32 



1 



I ^ n h .oh , h 

I - log"^ h log h log — 

327r^ V 3 32 32 32 



1 + 



41 



5127r3 



C(3)(B.167) 



We finally get 



(Pk 



lat (27r) 



A,ik)Aoik) 
' rC(3) 



d^k 1 



487r' 



Bvr^ Jfc2<32 (27r)2 k'^ 
647r 47r 6 4 



log TTT: - 1 log T^TT 



32 



32 



(B.168) 



We have checked numerically ( |B.163|) , ( |B.164|) and ( |B.168|) . This provides also a 
check of the integrals ( |7rT2q ) and ( |A.121| ). 

Using ( [B.146I) , ( [B.147|) and (|B.148| ) we finally get for the coefficients 02, ^2 and C3 



^ - l-^ + ^(K + 4V2-2V3)+0(l/iV) 



N 

_C3_ 



L + ^^K + W2-2Vs)+0{l/N) 

1 K Vi 

-^^(3) - - + ^ + 1^ + 0(l/iV) 
9o7r'^ 16 12 8 



(B.169) 
(B.170) 
(B.171) 



From these expressions, using the renormalization group relations ( ^.331 ) / ( p.41| ), we 
can compute the large- contribution to Wg"** and These expressions agree with 

(PD/dH). 

From the general expressions ( |B.140[ ) and ( |B.141D it is also possible to compute 
the value of higher-loops coefficients. Defining 



N' 



n-l 



(B.172) 
(B.173) 



we report their values in Table@. We can also compute, in view of the possibility 
of using improved expansions, the perturbative expansion of the isovector energy 



-et 




N- 1 




4(3 


~ 4iV 


N-l 


N 


4/3 


32(3 
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d^k 



1 iV" 



(B.174) 



n=3 



n 


an 


K 




fn 


1 


- 0.01413 


0.06259 


- 0.005993 


0.13650 


2 


- 0.01289 


- 0.01079 


- 0.001568 


0.02651 


3 


- 0.00404 


- 0.00401 


- 0.000530 


0.00929 


4 


- 0.00145 


- 0.00146 


- 0.000221 


0.00411 


5 


- 0.00063 


- 0.00063 


- 0.000109 


0.00212 


6 


- 0.00031 


- 0.00030 


- 0.000063 


0.00127 



Table 3: Values of the perturbative coefficients in the large- limit. 



We can thus rewrite for the correlation length 
where pE = {N - 1)/(4(1 - Ev)) and 
The coefficients and /„ are reported in Table (||). 
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